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Abstract 

The proper time path integral representation is derived explicitly for Green's func- 
tions in QCD. After an introductory analysis of perturbative properties, the total 
gluonic field is separated in a rigorous way into a nonperturbative background and 
valence gluon part. For nonperturbative contributions the background perturbation 
theory is used systematically, yielding two types of expansions,illustrated by direct 
physical applications. As an application, we discuss the collinear singularities in 
the Feynman-Schwinger representation formalism. Moreover, the generalization to 
nonzero temperature is made and expressions for partition functions in perturbation 
theory and nonperturbative background are explicitly written down. 
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1 Introduction 

The present stage of development of field theory in general and of QCD in 
particular requires the exploitation of nonperturbative methods in addition 
to summing up perturbative series. This calls for specific methods where the 
dependence on vacuum fields can be made simple and explicit. A good ex- 
ample is provided by the so-called Fock-Feynman-Schwinger representation 
(FSR) based on the Fock-Schwinger proper time and Feynman path integral 
formalism [1,2]. For QED asymptotic estimates the FSR was exploited in Ref. 
[3]. Later on this formalism was used in Ref. [4] for QCD and rederived in the 
framework of the stochastic background method in Ref. [5] . (For a review see 
also Ref. [6].) 

More recently some modification of the FSR was suggested in Ref. [7]. The 
one-loop perturbative amplitudes are especially convenient for the FSR. These 
amplitudes were extensively studied in Ref. [8] and a convenient connection 
to the string formalism was found in Ref. [9]. Important practical applica- 
tions especially for the effective action in QED and QCD are contained in 
Ref. [10]. Moreover, the first extension of FSR to nonzero temperature field 
theory was done in Refs. [11,12]. This forms the basis of a systematic study 
of the role of nonperturbative (NP) configurations in the temperature phase 
transitions [11, 13]. 

One of the most important advantages of the FSR is that it allows to reduce 
physical amplitudes to weighted integrals of averaged Wilson loops. Thus the 
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fields (both perturbative and NP) enter only through Wilson loops. For the 
latter case one can apply the cluster expansion method[14], which allows to 
sum up a series of approximations directly in the exponent. As a result we 
can avoid the summation of Feynman diagrams to get the asymptotics of form 
factors[15]. The role of FSR in the treatment of NP effects is more crucial. 
In this case one can develop a powerful method of background perturbation 
theory[16] treating the NP fields as a background[17]. 

In the present paper some of these problems will be treated systematically 
and in full detail, yielding a overall picture of the role of FSR in QCD. We 
will in particular focus on the relationship between the standard perturbative 
expansion and the FSR based expansion which clarify the important role of 
nonperturbative configurations in the vacuum. The previous publication of 
the authors on FSR in Ref. [18] was devoted to QED and cp 3 , cp 4 theories, and 
many basic formulas of FSR are already contained there. The later develop- 
ment of FSR in the framework of the field theory can be found in Ref. [19]. It 
has in particular been used successfully to reconstruct exact solutions of the 
one and two-particle Greens' function for <p 3 theory and scalar QED in the 
quenched approximation [20-25]. A review of the FSR applications to pertur- 
bation theory in QCD and a discussion of the connection between FSR and 
world-line formalism of Refs. [8-10] can be found in Ref. [26]. 

In the next section we describe how to derive the FSR formalism for the case of 
QCD. In section 3 we discuss the relationship between the usual perturbative 
expansion and the FSR. Section 4 deals with the study of two ways to deter- 
mine the Green's function, depending on the physical situation. One consists 
of an expansion in the perturbative fields and the other one is treating the 
nonperturbative fields as a correction. As applications we address in sections 
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5 and 6 the problems of collinear singularities and the finite temperature field 
theory in the FSR formalism, while some concluding remarks are made in the 
last section. 



2 General form of FSR in QCD 

Let us consider a scalar particle if (e.g. a Higgs boson) interacting with a 
nonabelian vector potential A, where the Euclidean Lagrangian is given by 



L v = ^\D^\ 2 + im 2 M 2 = \\{d^~ igAM 2 + t;™ 2 M 2 , (1) 



Using the Fock-Schwinger proper time representation the two-point Green's 
function of tp can be written in the quenched approximation as 

oo 

G(x,y) = (m 2 - D%1 = (x\P J dse^-^ly) . (2) 



To obtain the FSR for G a second step is needed. As in Ref. [1] the matrix 
element in Eq. (2) can be rewritten in the form of a path integral, treating 
s as the ordering parameter. Note the difference of the integral (2) from the 
case of the Abelian QED treated in Refs. [1,3,8]: A^ in our case is a matrix 
operator A^(x) = A a ^(x)T a . It does not commute for different x. Hence the 
ordering operator P in Eq. (2). The precise meaning of P becomes more clear 
in the final form of a path integral 

OO / X \ 

G(x,y) = Jds(Dz) xy e~ K Pexp iigjA^z)dzA , (3) 
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where K = m 2 s + \ Jq dr(-^) 2 . In Eq. (3) the functional integration measure 
can be written as 

(Dz) xy * lim ft IP^ f^M^^-y)) ( 4 ) 
v >xy n^oo A = y (4ne) 2 J (2vr) 4 KJ 



with Ne = s. The last integral in Eq. (4) ensures that the path z M (r), < r < 
s, starts at z^(0) = and ends at z^s) = x^. The form of Eq. (3) is the same 
as in the case of QED except for the ordering operator P which provides a 
precise meaning to the integral of the noncommuting matrices (zi), A^ 2 (z 2 ) 
etc. In the case of QCD the forms (3) and (4) were introduced in Refs. [4,5]. 

The FSR, corresponding to a description in terms of particle dynamics is 
equivalent to field theory, when all the vacuum polarisation contributions are 
also included [4, 5, 10], i.e. 



00 1 N rds r / r 

E Tn II /— (Dz t ) xx exp(-K) Pexp [ig A,(z)dz 

N = ^ • 1 = 1 J S * J \ y 

= J Dipexp (—Jd 4 xL v (x)^j . (5) 

Both sides are equal to vacuum-vacuum transition amplitude in the presence 
of the external nonabelian vector field and hence to each other. For practical 
calculations proper regularization of the above equation has to be done. The 
field A^ in Eq. (1) can be considered as a classical external field or as a 
quantum one. In the latter case the Green's functions (A.. A) induce nonlocal 
current-current interaction terms in the l.h.s. of Eq. (5). Such terms can also 
be generated by the presence of a (/9-field potential, in the r.h.s. of 

Eq. (5). 
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The advantage of the FSR in this case follows from the very clear space-time 
picture of the corresponding dynamics in terms of particle trajectories. This 
is especially important if the currents can be treated as classical or static (for 
example, in the heavy quark case). The mentioned remark on usefulness of 
the FSR (3) becomes clear when one considers the physical amplitude, e.g. 
the Green's function of the white state tr((p + (x)ip(x)) or its nonlocal version 
tr[ip + (x)<&(x, y) (p(y)], where $(x,y) - to be widely used in what follows - is 
the parallel transporter along some arbitrary contour C (x, y) 

$(x,y) = Pexp ^igjA^dz^j . (6) 



One has by standard rules 



G<p(x,y) = (tr \ip + (x)ip(x)\ tr \ip + (y)ip(y) 

oo oo 

= Jd Sl Jds 2 (Dz) xy (Dz') xy e- K - K ' (W) A + 



(7) 



o o 



where dots stand for the disconnected part, (G v (x, x)G<p(y, v))a- We have used 
the fact that the propagator for the charge-conjugated field (p + is proportional 
to &(x,y) = $(y,x). Therefore the ordering P must be inverted, &(x,y) = 
P exp(ig f% A fJi (z)dz^ l ). Thus all dependence on in G 9 is reduced to the 
Wilson loop average 



( W )a = (tr P c expigjA^z)dz^J 

C A 



(8) 



Here Pq is the ordering around the closed loop C passing through the points 
x and y, the loop being made of the paths z^(t), z' (t 1 ) and to be integrated 
over. 
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The FSR can also be used to describe the quark and gluon propagation. Similar 
to the QED case, the fermion (quark) Green's function in the presence of an 
Euclidean external gluonic field can be written as 

G q (x,y) = {^x)^(y)) q = (x\{m + by^y) 
= {x\(m - D)(m 2 - D 2 )- l \y) 

oo 

= (m-D) fds(Dz) xy e- K ^ a (x, y) , (9) 



o 



where <1> CT is the same as was introduced in Ref. [1] except for the ordering 
operators Pa, Pf 

$ a (x,y) = P A exp ^ig j A^dz^ P F exp J dra^F^ (10) 

with F/ju, = d^A v - d v A lx - ig[A^ A v ] and = ^(7//)V - while K and 

{Dz) xy are defined in Eqs. (3) and (4). Note that operators Pa, Pf in Eq. (10) 
preserve the proper ordering of matrices A^ and a^F^ respectively. Explicit 
examples are considered below. 

Finally we turn to the case of FSR for the valence gluon propagating in the 
background nonabelian field. Here we only quote the result for the gluon 
Green's function in the background Feynman gauge[5,17]. We have 

G^(x,y) = {x\{D\5, v - 2igF, v y 1 \y) (11) 

Proceeding in the same way as for quarks, we obtain the FSR for the gluon 
Green's function 

00 

G^(x,y) = J ds(Dz) xy e~ K ^^(x,y), (12) 


where we have defined 
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Now in the same way as is done above for scalars in Eq. (7), we may consider 
a Green's function, corresponding to the physical transition amplitude from a 
white state of q±, q\ to another white state consisting of qs, q^. It is given by 

G r q9 (x, y) = (G q (x, y)TG q (x, y)T - G q (x, x)FG q (y, y)T) A , (14) 

where T describes the vertex part for the interaction between the q, q pair in 
the meson. The first term on the r.h.s. of Eq. (14) can be reduced to the same 
form as in Eq. (7) but with the Wilson loop containing ordered insertions of 
the operators a^F^ (cf. Eq. (10)). 

3 Perturbation theory in the framework of FSR. Identities and 
partial summation 

In this section we discuss in detail how the usual results of perturbation theory 
follow from FSR. It is useful to establish such a general connection between 
the perturbation series (Feynman diagram technique) and FSR. At the same 
time the FSR presents a unique possibility to sum up Feynman diagrams in 
a very simple way, where the final result of the summation is written in an 
exponentiated way[15,17]. This method will be discussed in the next section. 

Consider the FSR for the quark Green's function. According to Eq. (9), the 
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2-nd order of perturbative expansion of Eq. (10) can be written as 

oo oo oo 

G q (x,y) = (m - D) jds Jdn Jdr 2 e~ K (Dz) xu d A u(Dz) uv d 4 v(Dz) vy 



x (igA^(u)u^ + ga^F^iy)) {igA u (v)v u + gcr X aFx a (v)) , (15) 

where we have used the identities 

(Dz) xy = (Dz) xu ^ l) d\(T 1 )(Dz) u(Tl)viT2) d 4 v(T 2 )(Dz) v(T2)y , (16) 

OO S T\ OO OO OO 

Jdsj dT 1 Jdr 2 f(s, n, r 2 ) = JdsJ dn Jdr 2 f(s + n + r 2 , n + r 2 , r 2 ). (17) 







We can also expand only in the color magnetic moment interaction (crF). This 
is useful when the spin-dependent interaction can be treated perturbatively, 
as it is in most cases for mesons and baryons (exclusions are Goldstone bosons 
and nucleons, where the spin interaction is very important and interconnected 
with chiral dynamics). In this case we obtain to the second order in (crF) 



oo oo oo 



G^{x,y) = %{m - b)JdsJdr l JdT 2 e- m " (s+Tl+T2) - K "- Kl - K2 



(Dz) xu ®(x, u)g(aF(u))d A u(Dz) uv $(u, v)g(aF(v))d A v(Dz] 



vy 



In another way it can be written as 



Gf\x,y) = i(m - D)(m 2 q - D% l u d A u g(aF(u))(m 2 q - Dj)" 1 * 

xg(aF(v))(ml-Dl)- y \ (19) 

Here (m 2 — D 2 )' 1 is the Green's function of a scalar quark in the external 
gluonic field A^. This type of expansion is useful also for the study of small- 
x behavior of static potential, since the correlator (aF(u)aF(v)) plays an 
important role there. 
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However, in establishing the general connection between perturbative expan- 
sion for Green's functions in FSR and expansions of exponential $ CT in Eq. (10), 
one encounters a technical difficulty since the coupling constant g enters in 
three different ways in FSR: 

1. in the factor (m — D) in front of the integral in Eq. (9) 

2. in the parallel transporter (the first exponential in Eq. (10)) 

3. in the exponential of g(crF). 

Therefore it is useful to compare the two expansions in the operator form 

(m + by 1 = (m + d- igA)' 1 = (m + By 1 + (m + <9) _1 i#A(m + d)' 1 
+(m + dy 1 igA(m + dy 1 igA(m + d)~ l + ... (20) 

and the FSR 

oo 

(m + by 1 = (m - D)(m 2 - d 2 )' 1 ]T(5(m 2 - d 2 y l ) n , (21) 

n=o 

where we have introduced 

5 = -ig(Ad + dA)-g 2 A 2 = D 2 -d 2 . (22) 

To see how the expansion (21) works, using D — d — igA, Eq. (21) becomes 

oo 

(m + by 1 = [(m + By 1 + %gk{m 2 - d 2 )' 1 ] £ [$(m 2 - d 2 )" 1 ]" 

n=0 

Separating out the first term we may rewrite this as 
{m + by 1 = (m + dy 1 

oo 

+ (m + d)-HgA(m - D)(m 2 - d 2 )' 1 £ [5(m 2 - d 2 )" 1 ]^) 

n=0 
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The last three factors in Eq. (23) are the same as occurring in Eq. (21). As 
a consequence the formal iteration of the resulting equation for the Greens' 
function reproduces the same series as in Eq. (20), showing the equivalence of 
the two expansions. 

It is important to note that each term in the expansion in powers of 5, after 
transforming the operator form of Eq. (21) into the integral form of FSR, 
becomes an expansion of the exponential $ CT in Eq. (10) in powers of g. The 
second order term of this expansion was written down before in Eq. (15). 



It is our purpose now to establish the connection between the expansion (21), 
(23) and the expansion of the exponential § a in Eq. (10) in the quark propa- 
gator (9). We can start with term linear in A and write (for the Abelian case 
see Appendix B of Ref. [18]) 



=ig J G^\x,z{T 1 ))d*z i -^A,{T 1 )Gf\z{T 1 ),y) ) 



(24) 



where the notation is clear from the general representation of G q , given by 
Eq. (9) 



N 



G g (x,y)=Jdse-^ll^ 2 e W 



f e(n) 



(25) 



with f (n) = z(n) - z(n - 1), A^n) = \[A^{z{n)) + A^z(n - 1))] and 



N 



= Pexp{ig £ Mn)ZM + 9 £ <V^(*( W M- (26) 

71=1 71=1 



Representing in Eq. (24) as |(^(L) +^(R)), where ^(L) refers to the 
integral over £ M in G^ to the left of £ M in Eq. (24) and ^(R) to the integral 
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in G^> standing to the right of £ M , we obtain 

= -,A e v. = 2iR ^, (27) 

Thus Eq. (24) in momentum space becomes 

= -9Gf{q){ q \p,A, + A, P ,\ q ')Gf{q') (28) 

In a similar way the second order term from the coinciding arguments yields 
Gf{coinc) = -g 2 J Gf\x, z)A%z)d i zGf (z,y). (29) 

Finally, the first order expansion of the term a^F^ in Eq. (10) yields the 
remaining missing component of the combination 5, Eq. (22), which can be 
rewritten as 

5 = -ig{A^ + d^) - g 2 Al + ga^F^. (30) 

Hence the second term in the expansion (21) 

(m + by 1 = (m- D)(m 2 - d 2 )' 1 

+ (m - D)(m 2 - d^dim 2 - d 2 )" 1 + ... (31) 

is exactly reproduced by the expansion of the FSR (9), where in the first 
exponential <3> CT in Eq. (9) one keeps terms of the first and second order, 0(gA^) 
and O^gA^) 2 ), while in the second exponential one keeps only the first order 
term O^ga^F^). It is easy to see that this rule can be generalized to higher 
orders of the expansion in 5 in Eq. (21) as well. 
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4 Perturbative vs nonperturbative: two types of expansion 



As was discussed in Section 2, gluons can also be considered in FSR. To make 
this statement explicit and to prove Eqs. (11-13) written for the gluon Green's 
functions, we can use the background perturbation theory [16]. As in Ref. [17] 
we combine the perturbative field a M and NP degrees of freedom in one 
gluonic field A^, namely 



At this point we must distinguish two opposite physical situations, which 
require different types of expansions. Consider first systems of small size, e.g. 
heavy quarkonia, which are mostly governed by the color Coulomb interaction 
and have a radius of the order (ma s (m)) _1 , where m is the quark mass. For the 
ground state bottomonium this radius is around 0.2 fm and for charmonium 



In this case we have the first type of expansion: at the zeroth order all gluon 
exhanges are taken into account (In practice the Coulomb contribution and 
the few first radiative corrections), while in first order one treats the nonper- 
turbative contribution as a correction. This expansion is considered in detail 
at the end of this Section. 

The second type of expansion takes into account NP interaction fully through 
NP vacuum correlators already in the zeroth order -this is the NP background 




(32) 



Under gauge transformations A^ transforms as 




(33) 



0.4 fm. 
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and in the next orders the usual background perturbation theory[16,17] is 
developed with necessary modifications. 



4-1 Expansion in perturbative fields 

We start this Section with the second type of expansion with some modifica- 
tions due to the independent integral over the background field, as in the 't 
Hooft's identity [17]. It is convenient to impose on the background gauge 
condition[16] 



In this case a ghost field has to be introduced. Defining = d\ ■ 5 ca + 



+ gf abc B b ^ = 0. 



(34) 



g f cba B b x = D\, We can write the resulting partition function as 




(35) 



where 



/fiC a r 1 

Da det(jj) exp J d 4 x[L + L(a) - — (G a ) 2 + J*a*]. (36) 



In Eq. (36) we have 




(37) 



and 



L(a) = Li(a) + L 2 (a) + L int (a) 



(38) 



with 
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L 1 {a) = alD c ;{B)F; u 

L 2 (a) = +-a v {D\8 tu/ - D^D V + igF^)a^ = 

= l -al[D™Df8» v - D™Df - g f^F^af, , 
L mt = - l -g {D,(B)a u - D u {B)a^) a f^a^ - ±g* 'f abc ayj ae - f a^(39) 

G a in Eq. (36) is the background gauge condition 

G a = d,al + gf abc BX = {D,a,f. (40) 
The ghost vertex is obtained from f£ = (D ll (B)D /Jl (B + a)) ab [16] to be 

L gh ost = -9+{D,{B)D,{B + a)) a %, (41) 

9 being the ghost field. The linear part of the Lagrangian L\ disappears if 
Bfj, satisfies the classical equations of motion. Here we do not impose this 
condition on B^. However, it was shown in Ref. [17] that L\ gives no important 
contribution 



We now can identify the propagator of a M from the quadratic terms in the 



Lagrangian L 2 (a) — ^(G a ) 2 . We get 



G% = [£>a<V " AA + igF» v + -D„b»]-J. (42) 

It will be convenient sometimes to choose £ = 1 and end up with the well- 
known form of the propagator in - what one would call - the background 
Feynman gauge 

G% = (D 2 X ■ V - 2igF llv )- 1 (43) 

This is exactly the form of gluon propagator used in Eq. (11). Integration over 
the ghost and gluon degrees of freedom in Eq. (36) yields 
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Z(J,B) = const{detW{B))r^ 2 \det{-D,{B)D,{B + a)} a=i _ 

oo c ( _S_\l I 

x{l + E ,, sj) exp(~JGJ) j v=D ^ )FM b) , (44) 
1=1 L - z 

where S int is the action corresponding to L(a) and G is defined in Eq. (43). 



Let us mention the convenient gauge prescription for gauge transformations 
of the fields a M , B^. Under the gauge transformations the fields transform as 

a M -> U+a^U, (45) 

B^U + (B^+ l -d^)U. (46) 

All the terms in Eq. (36), including the gauge fixing one \(G a ) 2 are gauge 
invariant. That was actually one of the aims put forward by 't Hooft in Ref. 
[16]. It has important consequences: 

(i) Any amplitude in the perturbative expansion in ga^ of Eqs. (36) and (44) 
corresponding to a generalized Feynman diagram, is separately gauge invariant 
(for colorless initial and final states of course). 

(ii) Due to gauge invariance of all terms, the renormalization is specifically 
simple in the background field formalism [16], since the counterterms enter 
only in gauge-invariant combinations, e.g. F 2 . The Z-factors Z g and Zb are 
connected: Z g Z B ^ = 1. 

As a consequence, the quantities like gB^, gF^ u (B) are renormalization-group 
(RG) invariant. Consequently all background field correlators are also RG 
invariant and they can be considered on the same footing as the external 
momenta in the amplitudes. This leads to a new form of solutions of RG 
equations, where a s = « s (Mb), Mb being the (difference) of the hybrid ex- 
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citations, typically Mb ~ lGeV. As a result a new phenomenon appears, 
freezing or saturation of a s at large Euclidean distances. For more discussion 
see Ref. [17] and recent explicit extraction of the freezing a s from the spectra 
of heavy quarkonia[31]. In the rest of this subsection we demonstrate how the 
background perturbation series (44) works for the meson Green's function. To 
this end we use Eq. (14) and consider the flavour nonsinglet case to disregard 
the second term in Eq. (14). 

Let us start with the meson Green's function and use the FSR for both quark 
and antiquark. 

oo oo 

G M (x,y) = (trT (f \m-D) J ds J dse~ K ~ R 



(Dz) xy (Dz) xy T^(m - D)W F ) (47) 
Here the barred symbols refer to the antiquark and 

s s 

W F = P A P F esxp(ig j dz^eMd J draft F^) exp(-^ j a$ F^df) .(48) 

o o 



In Eq. (47) integrations over proper times s, s and r, r occur, which also play 
the role of an ordering parameter along the trajectory, z^ = z^t), z^ = z^ij). 

It is convenient to go over to the actual time t = z± of the quark (or antiquark), 
defining the new quantity /x(t), which will play a very important role in what 
follows 

2/i(t) = ^, t = z 4 (r). (49) 
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For each quark (or antiquark and gluon) we can rewrite the path integral (47) 
as (see Refs. [38,42] for details) 

oo 

J ds(D 4 z) xy ... = const J D/i(t)(D 3 z) xy ... (50) 
o 

where (D 3 z) xy has the same form as in Eq. (4) but with all 4- vectors replaced 
by 3-vectors. The path integral D/i(t) is supplied with the proper integration 
measure, which is derived from the free motion Lagrangian. 

In general /x(t) can be a strongly oscillating function of t due to the Zitter- 
bewegung. In what follows we shall use the stationary point method for the 
evaluation of the integral over D/i(t), with the extremal f^oit) playing the role 
of an effective or constituent quark mass. We shall see that in all cases, where 
spin terms can be considered as a small perturbation, i.e. for the majority 
of mesons, /i is positive and rather large even for vanishing quark current 
masses m, rh, and the role of the Zitterbewegung is small (less than 10% from 
the comparison to the light-cone Hamiltonian eigenvalues, see Refs. [38,39] for 
details). 

Now the kinetic terms can be rewritten using Eq. (49) as 



T 



where T = rr 4 — y 4 . In the spin-dependent factors the corresponding changes 



are 



J dra^F^ = J ^—a^F^(z(t)). (52) 



Mi) 
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In what follows in this section we may systematically do a perturbation ex- 
pansion of the spin terms. They contribute to the total mass corrections of the 
order of 10-15% for lowest mass mesons, while they are much smaller for the 
high excited states. This perturbative approach fails however for pions (and 
kaons) where the chiral degrees of freedom should be taken into account. In 
this case another equation should be considered [40, 41]. 

Therefore as a starting approximation we may use the Green's functions of 
mesons made of spinless quarks. This amounts to neglecting in Eqs. (47,48) 
the terms (m — D),(fh — D) and o^F^. As a result, we have 

G { $(x,y) = const J Dfi(t)Dji(t)(D 3 z) xy (D 3 z) xy e- K ~ R (W). (53) 

The Wilson loop in Eq. (53) contains both perturbative and NP fields. It can 
be expanded as 

W(B + a) = W(B) 

oo 

+ J2(i9) n W {n \B ]X (l)...x(n))a^dx^(l)...dx^(n)- (54) 

n=l 

After averaging over a M , £> M we obtain, keeping the lowest correction term: 
(W(B + a)) B , a = (W(B)} B -g 2 (W {2 \B;x,y))dxdy + (55) 

where the second term in Eq. (55) can be written as 

- g 2 W {2) dxdy = -g 2 J $ a/3 (x, y, B) 

xt 5 a a CGf u (x, y, B)& s (y, x, B)dx,dy u , (56) 

G"^(x,y, B) being the gluon propagator in the background field (11). 
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We can easily see that W^> contains 3 pieces, two of them are perturbative self- 
energy quark terms. The third one, assuming that x and y refer to the quark 
and antiquark trajectory respectively, is the color Coulomb term, modified by 
the confining background. As argued in Ref. [17], this term represents the gluon 
propagating inside the world-sheet of the string between q and q. When the 
time T is large, the long film of this world-sheet does not influence the motion 
of the gluon, reducing it to the free OGE term. Hence (W^) factorizes into 
the film term ((W(B)) B ) and gluon propagator dx^dy^G yielding finally 
the color Coulomb term in the potential. (This is however only true for the 
lowest order term and only at large distances \x — y\ < T g ). Otherwise 
perturbative-nonperturbative interference comes into play [28]. 

In what follows we restrict our attention to the first term, (W(B)) b- Our next 
approximation is the neglect of perturbative exchanges in (W) (they will be 
restored in the final expression for Hamiltonian) . This yields for large Wilson 
loops, i.e. R,T > T g , 

(W) B = const exp(-aS min ) (57) 

where S m i n is the minimal area inside the given trajectories z(t),z(t) of the 
quark and antiquark, 



Here a point w on the surface is parameterized by = flz^t) + (1 — P)z^(t). 

The Nambu-Goto form of S min cannot be quantized due to the square root. 
To get rid of the square root we may use the auxiliary field approach [36] with 
functions v(/3,t) and rj((3,t) as is usually done in string theories. As a result 



T 



1 




(58) 



o 



o 
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the total Euclidean action becomes [42] 



A = K + K + aS mm = + g) + ^ 

, A*W -2 i -2 i / (T \2 2 o / • \ i 2 2i 



-r 2 + + (-) V - 2r?(wr) + ?? 2 r 2 ]}. (59) 
z, 2 z/ 

Here //+ = // + //, /j = ^gL, ^ = ^JfS r< = ^ - ^. Performing the 
Gaussian integrations over i? M and we arrive in the standard way at the 
Hamiltonian (we take m = fh for simplicity) 



L- r 2 a 1 !- 1 } dp , /• i/ 



—Jm + J2 d0 ' (60) 



^ + 2/„ 1 (/J-i)M/?M/3 2 J "(f) „ 
where p 2 . = (pr) 2 /r 2 and L is the angular momentum, L — (r x p). 

A reasonable approximation to the integrations over /i and z/ is to replace 
them by their corresponding extremum values[42]. For these values the terms 
fi(t) and z/(/5) have a simple physical meaning. E.g. when a = and L = 0, 
we find from Eq. (60) 



/J = 2 v /p 2 + m 2 , /i = yW™ 2 , (61) 

so that /i corresponds to the energy of the quark. Similarly in the limiting 
case L — > oo the extremum over z/(/3) yields 

= if 2 = 2iray/L(L + l). (62) 

V 1 - V(/3- 2) 

Hence z/ is the energy density along the string with f3 playing the role of the 
coordinate along the string. 
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4-2 Nonperturbative fields as a correction 



Now we turn to the first type of expansion mentioned above, i.e. when the NP 
contribution is considered to be a (small) correction to a basically perturbative 
result. As an example let us consider the spectrum of heavy quarkonia. We 
can calculate the NP shift of the Coulombic levels of heavy quark-ant iquark 
(qq) system, following Ref. [30]. When the quark mass m is large, the spatial 
and temporal extensions of the n-th bound state are 

n - n 2 , . 

, t n ~ -. (63) 



ma s ma z 



For low n ~ 1 these may be small enough to disregard the NP interaction in 
first approximation. So for the spin-averaged spectrum we can write 

M(n, I) = 2m{l- ^ + 0(a 3 s ) + A NP } (64) 

where Atvp is the expected nonperturbative correction, which should be small 
for states of small spatial extension. This conclusion can be drawn from the 
lattice (and phenomenological) parameterization of the static qq potential 

Act (r) 

V(r) = h or + const (65) 

3r 

Using the empirical values found for a = 0.2GeV 2 and a s (r) ~ 0.3 (at r m 
0.2/m) we may deduce that the first term on the l.h.s. of Eq. (65) is at 
r 0.3 fm comparable in magnitude to the second term, the NP contribution. 
Hence this suggests that the states with a radius r 0.3/m are mainly 
governed by the (color) Coulomb dynamics, while those with r ^> 0.3/m are 
mostly NP states. So we may expect e.g. the n — 1 bottomonium state to be 
largely Coulombic. 
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Let us now consider the general path-integral formalism for the qq system 
interacting via perturbative gluon exchanges and nonperturbative correlators. 
We start with the quark Green's function in the FSR form (cf. Eqs. (9,10)) 



S(x,y) = i(m- D) J ds(Dz) xy e~ K § a (x,y), 



(66) 



where 



K = m z s + - J zidr 
o 

and contains spin insertions into the parallel transporter 

X s 

$*(x,y) = P A PFexp[ig J A^dz fl + g J dra^F^z^))}. 



(67) 



Double ordering in and F^ is implied by the operators Pa, Pf- We have 
also introduced the 4x4 matrix in Dirac space 



P>i Ei 



(68) 



Neglecting spins, we have instead of Eq. (67) 

X 

$ a (x,y) -> $(x,y) = P A exp(ig J A^dz^. 



(69) 



In terms of the single quark Green's functions (66) and initial and final state 
matrices T^Tj ( such that g(x)Tj$(x, x)q(x) is the final qq state) the total 
relativistic gauge-invariant qq Green's function in the quenched approximation 
is similar to Eq. (14) 

G(x, x; y,y) = < trfifS^x, y)Ti<S>(y, y)S 2 (y, x)$(x, x))\ > 

- < tr(T f S 1 (x,xMx,x))tr(T t S 2 (y,y)^(y,y)) > . (70) 
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The angular brackets in Eq. (70) imply averaging over the gluonic field 
and the trace is taken over Dirac space. 

Since we are interested in this case primarily in heavy quarkonia, it is reason- 
able to do a systematic nonrelativistic approximation. To this end we introduce 
as in Refs. [34,42] the real evolution (time) parameter t instead of the proper 
time r in K, (f in K) and the dynamical mass parameters /i, p, as in Eq. (49) 

£ = = 2*; j^-fa^f. (7!) 



Here we have denoted 

T=]^{x A + x A ). (72) 

The nonrelativistic approximation is obtained, when we write for z±(t) , z±(t) 

z 4 (t) = t + c(t); z 4 (t)=t + C(t) (73) 



and expands in the fluctuations £, (, which are 0(^=). Note that the integra- 
tion in dsids 2 goes over into d/iidfj, 2 - Physically the expansion (73) means that 
we neglect trajectories with backtracking of z±, z^, i.e. dropping the so-called 
Z graphs. We can persuade ourself that the insertion of Eq. (73) into K, K 
allows to determine /xi, /i 2 from the extremum in K, K. We get 

^=7711 + 0(1/7711), n 2 = m 2 + 0(l/m 2 ). (74) 



We can further make a systematic expansion in powers of l/mj[34]). At least 
to lowest orders in 1/rrii this procedure is equivalent to the standard (gauge- 
noninvariant) nonrelativistic expansion [41]. 
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Let us keep the leading term of this expansion 

G(xx,yy) = 4m 1 m 2 e- (mi+m2)T J D 3 zD 3 ze~ Kl ' K2 < W(C) >, 



(75) 



where K x = ?f J T zf(t)dt, K 2 = ^ z^dt. Furthermore, < W(C) > is 
the Wilson loop operator with a closed contour C comprising the q and q 
paths, and the initial and final state parallel transporters <3>(x, x) and $(y, y). 

The representation (75) will be our main object of study in the remaining part 
of this Section. For the heavy qq system the perturbative interaction contains 
an expansion in powers of ^f-(v being the velocity in the cm. system). This 
should be kept entirely, while the nonperturbative interaction can be treated 
up to the lowest order approximation. For the total gluonic field we may 



where is the NP background, while a M is the perturbative fluctuation. 

It is convenient in this first type of expansion to split the gauge transformation 

as 



write 




B u + a 



(76) 



B. 



V + B,V, a 



V + (a, + -d,)V, 



(77) 



so that the parallel transporter 



X 




(78) 



y 



transforms as 



&(a;x,y) -> V + (x)$(a; x, y)V{y). 



(79) 
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Using Eq. (75) we may now determine the effects of the NP contribution as 
a correction. The Wilson loop average in Eq. (75) can be written, using Eq. 
(76), as 



< W (C) >=< trPexp(ig J A^dz^) >=< trPexp(ig J a^dz^) > 

c c 



(ig) 



j-Jdz^J dz'^ < trP<S>(a;z,z')B u (z')$(a; z',z)B^{z) > + ... 
c c 

W + W 2 + (80) 



where we have omitted the term linear in since it vanishes when averaged 
over the field B^. The dots imply terms of higher power in B^. The contour 
and points z,z' are schematically shown in Fig. 1. 




Fig. 1. The contour C, characterized by the quark trajectories z and z' in the Wilson 
loop with possible ladder-type gluon exchanges between the quarks. 

Let us discuss the first term on the r.h.s. of Eq. (80). It is the Wilson loop 
average of the usual perturbative fields, discussed extensively in Ref. [27]. We 
can use for Wq the cluster expansion to obtain 

W = Zexp(^ 2 + y?4 + y? 6 + ■■■), (81) 
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V>2 = 




// 




iv c 2 -i 

2iV c 



(82) 



where regularization is implied in the integral (f 2 to be absorbed into the Z 
factor. Note that (p 2 contains all ladder-type exchanges. In addition also the 
"Abelian-crossed" diagrams — those where times of the vertices can not be 
ordered while color generators tf are always kept in the same order, as in the 
ladder diagrams. Therefore all crossed diagrams (minus "Abelian-crossed") 
are contained in <^ 4 and contribute 0(l/iV c ) as compared with ladder ones 
(cf the discussion in Ref. [27]). In addition ip4 contains "Mercedes-Benz di- 
agrams", again repeated infinitely many times. It is interesting to note that 
each term <p 2 , <f4 etc. in Eq. (81) sums up to an infinite series of diagrams. 

In particular exp(ip 2 ) contains all terms with powers of ^f-, as we shall see 
below. For heavy (and slow) quarks we can write: 



where f=z — z f . In this way we obtain a singlet one-gluon-exchange (OGE) 
potential, the effective time difference being At = \t — t'\ ~ |f|. In addition 
also the radiative corrections due to the transverse gluon exchange can be 
obtained in this way. For the qq mass <p 2 leads to a correction of order 0(a s ), 
being the Coulombic energy. It should however be noted that in the wave- 
function the Coulomb potential has to be kept to all orders because of its 
singular character. In fact we shall not expand exp{tp 2 ^), while this is done 
for the other contributions like the radiative corrections. 





o 



(83) 
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Turning to W 2 we may determine the leading (in N c ) set of diagrams. They 
consist of the diagrams, where the gluons propagate between the q and q lines 
with the same time coordinates, i.e. diagrams with Coulombic or instantaneous 
gluon exchanges. We have for the gluon propagator 

< a > )a>M >= Jt-vr m 



Here we have chosen for simplicity the Feynman gauge, since W and W 2 are 
gauge invariant. Expanding W 2 in Eq. (80) in powers of a M we have in view of 
Eq. (84) terms typically of the form 

tr (t bk t hk -\ ..t bl t bl t bk t ai ...t an t a Bp an ...t ai B b /.. . ) 

-> (C 2 ) k tr(t ai ...t an t a Bp an ...t ai B b /...) (85) 

Now due to the equality 

m = ~wf> m 

We obtain for all exchanges in the time interval between times of B^(z) and 
B u (z'), a factor (— ^r) instead of the factor C 2 for all the other exchanges. 

The correction term W 2 can be worked out explicitly. In particular, we now 
derive the lowest order corrections to the energy levels and wave-functions 
due to the NP field correlators. As shown in Fig. 2 we may divide the total 
time interval T into three parts 

(1) < t < w' 4 

(2) w' A < t < w A 

(3) w 4 < t < T 
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Fig. 2. Contribution of the Gaussian correlator to the qq Green's function. The points 
w and w' of the correlator are connected by the parallel transporter, shown by the 
solid line going through the point xq- This makes the correlator gauge invariant. 

where t is the cm. time. In K + K we may separate the cm. and relative 
coordinates 



m 1 z i + m 2 z i 

Ri = ■ ; nit) = Zi{t) - Zi(t), (87) 

m 1 + m 2 



so that 



K + K 



T MR 2 , m } . , 

~Y~ + ~2 J (t) ' 

o o 



I 



(88) 



withm=^^, M = m 1 +m 2 . 
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Separating out the trivial cm. motion, we have in the parts (1) and (3) the 
path integrals, representing actually the singlet Coulomb Green's function 

G { c\r{h),r{t 2 )-h - t 2 ) = J Dr{t)e^^ dt + C2aa £l (89) 

In the part (2) instead we have an octet Coulomb Green's function 

GP(r(w 4 ),r(w' 4 );w 4 -<)= J Dr(t)e 2 } < 2Nc J < lml . (90) 

As a result W 2 can be written as 

W^^ffdz,/ K < trB v B p > X dtV °< V * d \ (91) 



c c 



where 



V ° = C 2 ^, VS = -— — . (92) 



We can easily identify — V® and — Vj as a singlet and octet qq potential, con- 
sidered in Ref. [30]. For < trB u B^ > we can use the modified Fock-Schwinger 
gauge to obtain: 

z 

B^z) = J dw p a(w)F pil (w) (93) 

and 

Jdz^dzl < B^B^z ) >= J da^da'^ < F^(w)F uX (w') >, (94) 

where we have introduced a surface element da w = dz^dw p a(w). To make Eq. 
(94) fully gauge-invariant we can introduce in the integral in Eq. (93) factors 
$(rr ,u>) identically equal to unity in the Fock-Schwinger gauge. As a result 
we get 
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< trB u B^ > =< J da w da v \ 

xtr{$(x , w)F[ 1 p(w)'&(w, xq)^(^xq, w ^)F v \{w , xo)}(9§,) 

which is fully gauge invariant. 

Using Eq. (91) we find for the correction to the total qq Green's function 
G = + AG, 



q 2 



da^(w) J da f ,, l/ ,(w')d 3 r(w 4 )d 3 r(w 4 ) 

G^(r(T),r(w 4 );T- w 4 )Gf\r(w 4 ),r(w' 4 ),w 4 - w' 4 ) 

x < F^{w)F^{w') > G^(r(w' 4 ),r(0);w' 4 ), (96) 

where the integrals over da^ u , da^ v i are taken over the surface E^. It is 
convenient to identify E^ with the minimal surface inside the contour C, 
formed between the trajectories z(t) and z(t). Introducing the straight line 
between z(t) and z(t) 



w,(t, (5) = z,{t)(5 + z,(t)(l -p)=R^ + r (/3 (97) 

mi + m 2 



with 

w 4 = z 4 = z 4 = t 
we may write the surface element as 

da^viw) = {w'^w u — w pw'^dtdfl = a^dtdfl (98) 



with < = ^ 



~~dt 



. We also have in the cm. system 



t 1 (a mi \ 

o-a — Ti, ciij — e.iik ij k~r^\P ; )i 

%m mi + m 2 



(99) 
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while the Minkowskian angular momentum L is given by 
1 d 

Li = e ikl r k ■ - — . (100) 

l OTi 



In the nonrelativistic approximation we expand in powers of — , Hence 
can be neglected in lowest order and we are left with only a i4 , i.e. in Eq. (96) 
only the electric field correlators should be kept. The field correlators have the 
following representation in terms of the two Lorentz invariants D and D\ [29] 

g 2 tr < Ei(w)E k (w') >= 

^[S ik (D(w - w') + Dxiw - w') + hf^) + KhJ^l (101) 

where hi = Wt — w[. D and D 1 are normalized as 



£>(0) + DM =g 2 < trF%(0) >= jlrrf (102) 



G 2 is the standard definition of the gluonic condensate [32] 



G2 = ^ < F; V F; U >= O.OUGeV 4 . (103) 



Inserting Eq. (101) into Eq. (96) and neglecting the terms h;Jn k ~ 0(4a) we 
get 



AG = -—G { c\r(T),r) d 3 rG { c\r,r')d 3 r' 

r i dpdw 4 ! J i dpdw' A &(w - w')G^(r', r(0)), (104) 



where we have defined 



3D 

A(w - w') = D(w - w') + D^w - w') + h\-^. (105) 
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Using the spectral decomposition for G c 



G^ % \r,r',t) =< r\e- H ' mt \r' >= V ^' 8 \r)e- E -'%^ + (r'), (106) 



we can rewrite Eq. (96) for the matrix element of AG between singlet Coulomb 
wave functions 



_ F (i) T , , ._ 
e J_ T [ dp 4 dp ? 

24 J (2tt) 4 
< n r,e ra i+ m 2 / k >< k r'e m i+ m 2' \n> 



<n\AG\n>= — — T / gA(p)d^ 



(2tt)< 

(8) 



E — ( ;: — . ( lor ) 

fc=o,i,... E k E n- m 



where A(p) is the Fourier transform of Eq. (105). The set of states \k > in Eq. 
(107) with eigenvalues i?^ refer to the octet Hamiltonian piece in Eq. (90) 

Jf < 8 > = i + (108) 
2m 2iV c |r1 v ; 



The correlator A(x) depends on x as A(x) = f(jr), and decays exponentially 
at large |x|[42]. For what follows it is crucial to compare the two parameters, 
T„ and the Coulombic size of the n-th state of the qq system, R„ = — . In 

y mc20L s 

the Voloshin-Leutwyler case [34] it is assumed explicitly or implicitly that 

Case (i) T g > R n 

In the opposite case 

Case (ii) T g <C R n 
as we shall see completely different dynamics occurs. 
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Writing G + AG = const e -( E ^ +AE ^ T w const e~ E ™ )T (l — AE n T) we finally 
obtain in case (i) for AE n 



AP ir 2 G 2 < n\rj\k >< k\rj\n > 

= — 8 • ( 109 ) 



Results of the calculations [30] using the Voloshin-Leutwyler approximation 
(VLA) (Eq. (109) for charmonium and bottomonium) and using the standard 
value of G 2 [32] are given in Table I. One can see that a rough agreement 
exists only for the lowest bottomonium state. Consider now the opposite case, 
T g <C R n . Since lattice measurements yield T g ps 0.2/to, we may expect that 
this case is generally applicable to all bb and cc states. However in this case it is 
not enough to keep only the W2, but also sum up all NP terms, which amounts 
to the exponentiation of the NP contribution. The NP local potential appears 
in addition to the Coulomb term. This has been done in the framework of the 
local potential picture in Ref. [35], where the NP potential is expressed via 
correlators D(x) and D 1 (x). 

Results of calculations of Ref. [35] yield a very consistent picture both for 
levels and wave functions of bottomonium and quarkonium. To compare with 
VLA and our results here the results of Ref. [35] are listed in the middle 
column of the Table, demonstrating a much better agreement with experiment 
than results for the VLA. Note, that in Ref. [35] the NP interaction was not 
treated as a perturbation, but nonperturbatively by including the NP part in 
the potential. Consequently this explains why the predictions have improved. 
To summarize, because of the small T g ps 0.2/m, the potential picture is more 
adequate for quarkonia than the VLA formalism or QCD sum rules, including 
even the bottomonium case. 
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Table 1. 

The experimental values and the predicted splittings in MeV of various 
states in bottomium and charmonium in the Voloshin-Leutwyler 
approximation and Ref. [35]. 

splitting VLA [35] exp. 

2S-lS(bb) 479 554 558 

2S-2P(bb) 181 112 123 

3S-2S(bb) 4570 342 332 

2S-lS(cc) 9733 582 670 



5 IR and collinear singularities in FSR 

It is known that in QED some matrix elements and partial cross-sections dis- 
play singularities [43], which are of two general types; a) due to soft photon 
exchange (IR singularities) b) due to the collinear motion of an exchanged 
and emitted photon (collinear singularities). A similar situation exists in per- 
turbative QCD, where the same lowest order amplitudes contain both IR and 
collinear singularities, see Ref. [44]. 

Examples of Feynman graphs which produce both types of singularities are 
given in Fig. 3. These graphs refer to the process e + e~ — > qq and the total 
cross section for the sum of three graphs where an additional gluon can be 
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(a) (b) (c) 

Fig. 3. Feynman graphs to order a s with singularities, which are cancelled in the 
crossection. 

generated is finite to the order 0(a s ) and is equal to 

^ (1) = ^ (0) {i + ^^}- (HO) 

This is in line with the KLN theorem[45] derived and proved in the framework 
of QED. 




Fig. 4. The gauge-invariant 7 — 7 amplitude, where all sigularities cancel, while they 
are separately present in the parts of graphs obtained by cutting along the dash-dotted 
lines. 

We now argue that the situation is different in QCD when the nonperturbative 
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confining vacuum is taken into account. This in particular can be demonstrated 
in the framework of the FSR. We start with the QCD contribution to the 
photon self-energy part, which is shown graphically in Fig. 4. It has the form 

IWQ) = (Q,Qu - Q\ V )U(Q 2 ). (Ill) 



In particular, the 0(a s ) part of IT^ (Q 2 ), shown in Fig. 4b, yields an analytic 
function of Q 2 , with the imaginary part (discontinuity across the cut Q 2 > 0) 
given by the sum of the three graphs of Fig. 4a,b multiplied by the comple- 
mentary parts as shown in by the dash-dotted lines in Fig. 4b. It is clear, that 
both the whole function H(Q 2 ) and its total absorptive part is free from IR and 
collinear singularities, while each piece in the imaginary part, yielding partial 
crossections cr( ), cr(p,c) corresponding to the graphs Fig. 3a-c are IR divergent. 
At this point the difference between QED and QCD can be felt even on the 
purely perturbative level. Namely, in QED the process e + e~ — > e + e~ + ny 
cannot be associated with the imaginary part of photon self-energy li^Q 2 ), 
since photons can be emitted in any amount off the il 7 (Q 2 ) and hence should 
be summed up separately. 

This fact is formulated as a notion of a physical electron, containing a bare 
electron plus any amount of additional soft photons (see the Bloch-Nordsiek 
method discussed for example in Ref. [43]). In QCD the situation is different, 
since separate gluons cannot escape the internal space of li^Q 2 ) (cannot be 
emitted), except when they create (pairwise, triplewise etc.) massive glueballs, 
or else when they are accompained by the sea quark pairs forming hybrid 
states. 

To illustrate our ideas we shall use the FSR, introducing the background 
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confining field, and using the method of Ref . [47] . We define the photon vacuum 
polarization function H(q 2 ) as a correlator of electromagnetic currents for the 
process e + e~ — > hadrons in the usual way 

- i J d 4 xe 1 ^ < 0\T( Jfl (x) 3u (0)\0 >= {q^q v - 9flu q 2 )U(q 2 ), (112) 
where the imaginary part of II is related to the total hadronic ratio R as 

m = ^ r ^) = „t M (113) 

There are two usual approaches to calculate II (g 2 ). The first one is based on 
a purely perturbative expansion, which is now known to the order 0(0^) [50]. 
The second one is the OPE approach[32], which includes the NP contributions 
in the form of local condensates. For two light quarks of equal masses {m u = 
m d = m) it yields for U(Q 2 ) 

1 ,„ a s ., Q 2 6m 2 2m < qq > a s < FF > lS 

nw » = + - )ln iJ + 1? + + ~^Q^ + - (114) 

In what follows we shall include the NP fields as they enter into Green's 
functions, i.e. nonlocally. Moreover, we shall be mostly interested in the large 
distance behaviour, where the role of NP fields is important. To this end we 
first of all write the exact expression for II (Q 2 ) in the presence of the nonper- 
turbative background and formulate some general properties of the perturba- 
tive series. More specifically, the e.m. current correlator can be written in the 
form [17,47] 

n(Q 2 ) = 1 J e iQx d 4 x J DB J Dae- s ^ B+a ^ 

xtr(^G q (x, 0)7 M Gg(0, x))det(m + D(B + a)) (115) 
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Here G q is the single quark Green's function in the total field + a M , 

G q B+a \x,y) = <x \( m + d- ig{B + a))" 1 ^ > . (116) 

The background quark propagator is conveniently written using the FSR as 

OO X s 

Gf\x,y) = J dse- K DzPexp[tg J (B^dz^} exp[g J a^F^dr] (117) 
y 

With K = \j°Z 2 (T)dT. 

To simplify our analysis we take the limit N c — > oo end drop the det term 
in Eq. (115). We are then left with only planar diagrams containing gluon 
exchanges in the external background field. Moreover, all gluon lines 

G^ in the limit N c — > oo are replaced by double fundamental lines[51] and 
we are left only with diagrams, where the area S between the quark lines in 
H(Q 2 ) is divided into a number of pieces ASk, shown schematically in Fig. 5. 




Fig. 5. A generic Feynman diagram for the photon self- energy in high order of 
background perturbation theory at large N c . All areas ASk between double gluon 
lines are covered by the confining film, yielding the area law in Eq. (118). 

The infrared behaviour of a s (Q 2 ) at small Q 2 is connected to the limit of 
large areas of S. In this limit the product of all phase integrals = 
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exp{ig Bp.dzp) from all Green's functions can be averaged using the area 
law, i.e. we have (modulo spin insertions aF, which are unimportant for large 
distances[17]) 

< i^Mxi, yi ) > B = n™ =1 < w(As k ) > 

« exp(-(7 ]T AS k ) for N c -> oo. (118) 

fc=i 

This last factor serves as the IR regularizing factor in the Feynman integral, 
preventing any type of IR divergence. Using representations (115) and (117), 
we can formulate the following theorem. 

Theorem: 

Any term in the perturbative expansion of U(Q 2 ), Eq. (115), in powers of ga^ 
at large N c can be written as a configuration space Feynman diagram with 
an additional weight < W s (Ci) > for each closed contour. Here W s {Ci) is the 
Wilson loop with spin insertions, as in Eq. (117). Brackets denote averaging 
over background fields. 

This theorem is easily proved by expanding Eq. (115) in powers of ga^ and 
using double quark lines, Eq. (117), for gluon lines. Since at large N c we can 
replace adjoint color indices by doubled fundamental ones, then in each planar 
diagram the whole surface is divided into a set of closed fundamental contours, 
for which Eq. (118) holds true, again due to large N c . Thus to each contour is 
assigned the Wilson loop < W(Ci) >. The rest is the usual free propagators 
written in the FSR. 

Looking now at the large distance behaviour of the resulting configuration 
space planar Feynman diagram, we may derive from the above theorem the 
following corollary. 
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Corollary: 

Any planar diagram for II (Q 2 ) is convergent at large distances in the Euclidean 
space-time in the confining phase, when < W(Ci) >~ exp(-aSi). 

The proof is trivial, since the free planar diagram may diverge at large dis- 
tances at most logarithmically. The kernel < W(Ci) > makes all integrals 
convergent at large distances. At small distances (i.e. for small area, Si — > 0) 
the kernel < W(Ci) > behaves as 

< W(Ci) >~ exp( ) 

(see second reference in Ref. [29]). Hence the structure of small-distance per- 
turbative divergencies is the same for the planar diagram whether the NP 
background is present or not. Therefore the usual renormalization technique 
(e.g. the dimensional renormalization) is applicable. As a result the planar 
Feynman diagram contributions to a™U( n \Q 2 ) in the background is made fi- 
nite also at small distances. The consequence of this is that any renormalized 
term a™n( n )(Q 2 ) in the perturbative expansion of U(Q 2 ) is finite at all finite 
Euclidean Q 2 , including Q 2 = 0. 

Then it follows that we can choose the renormalization scheme for a s , which 
renders a s finite for all < Q 2 < oo and the Landau ghost pole will be absent. 
To make explicit this renormalization of a s , we can write the perturbative 
expansion of the function U(Q 2 ), Eq. (115), as 

n(Q 2 ) = n(°)(Q 2 ) + a s n«(Q 2 ) + ^n^Q 2 ) + ... (119) 

We now again use the large iV c approximation, in which case contains 
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only simple poles in Q 2 [51]: 

1 00 C 
n( ° )(g2)= 12^5 Q^f 

where the mass M n is an eigenvalue of the Hamiltonian 
quarks and background field B^, 

while the constant C n is connected to the eigenf unctions of [47] 

n ~ M n ' 

where 

2tt 2 7 E fc 3 (£* + m) 

00 

A 2 = 27r 2 (| dkuKk))- 1 ; E k = (k 2 + m 2 ) 1 ' 2 . 


In what follows we are mostly interested in the long-distance effective Hamil- 
tonian. It can be obtained from G q q for large distances, r » T 9 , where T g is 
the gluonic correlation length of the vacuum, T g ps 0.2/ra[33]: 

G g g(x,0) = I DB^tr^G^x^G^x)) =< x\e- H(0) \ x \\$ > . (123) 

At these distances we can neglect in Eq. (117) the quark spin insertions o^F^ 
and use the area law: 

< W c >^ exp(-aS min ), (124) 
where S m i n is the minimal area inside the loop C. 



(120) 



. It contains only 



(121) 



. We have 
(122) 
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Then the Hamiltonian in Eq. (123) is readily obtained by the method of Ref. 
[36]. In the cm. system for the orbital momentum I — it has the familiar 
form: 

# (0) = 2^Jp 2 + m 2 + ar + const, (125) 



where a constant appears due to the perimeter term in < Wc >• For / = 2 a 
small correction from the rotating string appears[42], which we neglect in first 
approximation. 

Now we can use the results of the quasiclassical analysis of H^[52], where the 
values of M n , C n have already been found. They can be represented as follows 
(n = nr + 1/2, n r = 0, 1, 2, I = 0, 2) 

Ml = 2vr<r(2n r + /) + M 2 , (126) 



where Mq is a weak function of the quantum numbers n r ,l separately, com- 
prising the constant term of Eq. (125). In what follows we shall put it equal 
to the p-meson mass, Mq ~ m 2 p . For C n one obtains quasiclassically[52] 

C n {l = 0) = ^Q)N c ml ml = 4na, 
C n (l = 2) = \Q)N c ml (127) 

Using the asymptotic expressions Eqs. (126-127) for M n , C n and starting with 
n = no, we can write 



n (o)^ 1 y Cn Q}N C Q 2 + M 2 + n ml 
W ' ~ 12vr2 n ^ Ml + Qi 12vr2 n m% ] 

+ divergent constant. (128) 
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Here we have used the equality 

^ 1 1 ,,Q 2 + M 2 + n m 2 

> — — — = Kipl 7, )+ diverqent constant (1 29) 

and ip(z) = 

In Eq. (128) we have separated the first n terms to treat them nonquasi- 
classically, while keeping for the other states with n > n the quasiclassical 
expressions (126-127). In what follows, however, we shall put n — 1 for sim- 
plicity. We shall show below that even in this case our results will reproduce 
e + e~ experimental data with good accuracy (see Ref . [47] for details) . 

Consider the asymptotics of Tl^(Q 2 ) at large Q 2 . Using the asymptotics of 



1 ^ B 



^)*-oc = lnz- — -Y, ' ( 13 °) 

k=l 



where B n are Bernoulli numbers, we obtain from Eq. (128) 



Q)N C Q 2 + M?, 



n «»» W ^-^l n ^ + 0(^). (1 31) 

We can easily see that this term coincides at Q 2 ^> Mg with the first term in 
the OPE (114) - the logarithmic one. Taking the imaginary part of Eq. (131) 
at Q 2 — > — s we find 

R(q 2 ) = 127r/mn (0) (-s) = N c Q 2 f , (132) 

i.e. it means that we have obtained for nJ ) the same result as for free quarks. 
This fact is the explicit manifestation of the quark-hadron duality. 

The analysis of more complicated planar graphs as in Fig. 5 can be carried 
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out as in Ref. [47], yielding in this way a new perturbative series with a s 
renormalized in the background fields and having therefore no Landau ghost 
poles. We refer to Refs. [17] and [47] for more details. 

Now instead we return to the diagrams in Fig. 3, which give the lowest order 
perturbative amplitudes associated with the production of 2 and 3 jets. The 
corresponding photon self-energy part is depicted in Fig. 4b. When the non- 
perturbative interaction is disregarded, even in the hadronization process, so 
that a gluon emitted with the moment k is directly associated with a gluon jet, 
then the singularities of the partial cross sections of Fig. 3 are transmitted into 
the singularities of the jet cross sections. These are cured by the introduction 
of the jet thickness, as in Sterman- Weinberg method[48] or introducing finite 
angular resolution rjo to distinguish 2-jet and 3-jet events (see Refs. [44], [53] 
for details). In what follows we show, that in the FSR with account of back- 
ground fields, all IR and collonear singularities disappear. Therefore the cross 
sections for 2-jet and 3-jet events remain finite. 

As was discussed above, in the leading l/N c approximation we have only 
planar graphs for T[(Q 2 ), describing "1-jet events", which actually create a 
constant behaviour of R(q 2 ) (apart from new opening thresholds), exactly 
reproducing the hadronic ratio, see Eq. (132). Speaking of 2-jet events we 
may actually consider the next approximation in 1/N C , since we need an extra 
quark loop in the Tl(Q 2 ). This can be easily derived from Eq. (115), where the 
determinant can be expanded in the FSR 

lndet(m + Z)) = ]- ln[det(m 2 - D 2 )] = \tr ln(m 2 - D 2 ), (133) 
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where we have used the symmetry property of the spectrum of D. Hence 



ds 



det(m + D)=expl--tr J — £(s)e~ sm ~ Dz xx W a (A, F) 



(134) 



d M 2s t s 



o or 



with £(s) a regularizing factor. For this we may take £(£) = lim^ r ^ }s 
we can use the Pauli-Villars form for Furthermore in Eq. (134) we have 



s 

W a (A,F) = P A P F expig J A^dz^ -expg J a^F^dr. 



(135) 



It is clear that det(m + D) allows an expansion in the number of quark loops, 
which is done by expanding the exponential in Eq. (134). Keeping only one 
quark loop for the 2-jet events, we obtain the graph, shown in Fig. 6 with 
an internal quark loop from the determinant. It is essential, that the whole 
region between the loops is covered by the NP correlators, creating a kind of 
"film"- the world surface of the string, with perturbative (i.e. generated by 
d/j) exchanges. 




Fig. 6. Photon self-energy graph corresponding to the 2-jet cross-section with one 
dynamical quark loop. 

It is clear that in this situation quarks are never on the mass shell (in contrast 
to the purely perturbative case) and therefore both IR and collinear singu- 
larities are absent. We can consider also the amplitude for the 3-jet event, 
with one gluon jet. Its perturbative amplitude corresponds to Fig. 3b, c. The 
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perturbative situation is discussed in detail in Ref. [53], and the 3-jet cross 
section for the process e + e~ — > qqg is given by 

I-ffe- = C F gl Xl l %l - , (136) 
a dxidx2 27T (1 — — X2) 

where the integration region is < x±, x 2 < 1, x 1 + x 2 > 1. The integral 
is divergent both due to collinear and IR effects, since 1 — x± — x 2 E g {l — 
cos 6 2g ) I \fs and 1 — x 2 = X\E g {\ — cos6>i 9 )/ \fs, where E g is gluon energy and 
6i g is angle between the gluon and i-th quark. To handle these divergencies 
one can use the so-called JADE algorithm [54], where the minimum invariant 
mass of a parton pair is larger than ys, i.e. min(pi + pj) 2 > ys. With this 
condition the energy region for 3-jet events looks like 

< x 1 ,x 2 < 1 - y, x l +x 2 >l+y. (137) 



The nonperturbative counterpart is obtained in two ways: a) the emitted gluon 
is accompanied by another gluon, forming together a two-gluon glueball (as 
was calculated in the framework of FSR in Ref. [49]), or b) a hybrid formation, 
when the emitted gluon is accompanied by a sea quark-antiquark pair. Both 
possibilities are depicted in Fig. 7. For the hybrid case we should expand 
the determinant term (the exponential in Eq. (134)) to the second power, 
producing in this way the two quark loops. 

It is clear in this case, that all particles, including the gluon, are off-shell 
and IR and collinear singularities are absent. Moreover, assuming as usual 
almost collinear hadronisation[53] we should replace the momenta of quarks 
and gluons by the corresponding momenta of hadrons. We can easily see that 
the factor is singular in case of qqg system becomes 2pi k+AM 2 ' wnere 
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(a) (b) 



Fig. 7. Graphs corresponding to the 3-jet cross-section with the gluon hadronized 
into a glueball or accompanied by a sea-quark pair forming a hybrid. 

AM 2 = M 2 - Mf + M 2 g and M g is the hybrid (glueball) mass, so that AM 2 
is in the GeV region. It effectively cuts off the singularity at small y cut , as 
is clearly seen in the experimental data[55]. The same reasoning applies for 
higher jet events. The only issue which exists is of experimental character. 
It amounts to the precise definition of the number of jets, i.e. to classify the 
hadrons between the several jets. 

Thus experiment as well as background perturbation theory do not exhibit 
collinear and IR singularities pertinent to the standard perturbation theory. 

6 FSR at nonzero temperature 

Within the framework of the FSR the problem of the various Green's func- 
tions at finite temperature can be studied[ll,12]. We first discuss the basic 
formalism for T > and then turn to the calculation of the gluon and quark 
Green's functions. 
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6.1 Basic equations 



We start with standard formulae of the background field formalism [16, 17] 
generalized to the case of nonzero temperature. We assume that the gluonic 
field can be split into the background field £> M and the quantum field a M 

A-jx — Bp + a M , (138) 

both satisfying the periodic boundary conditions 

B^z^Zi) = B^z* + nP,Zi); a M (z 4 , z { ) = a^(z 4 + n(3, z t ), (139) 

where n is an integer and f3 = 1/T. The partition function can be written as 
Z(V,T) = (Z(B)) B (140) 

with 

Z(B) =N J D0exp i^-j dT j d 3 xL tot (x,r)^j (141) 

and where denotes all set of fields a M , L tot is the same as L(a) defined 

in Eq. (38) and iV is a normalization constant. Furthermore, in Eq. (140) 
( ) B means some averaging over (nonperturbative) background fields B^. The 
precise form of this averaging is not needed for our purpose. 

Integration over the ghost and gluon degrees of freedom in Eq. (140) yields the 
same answer as Eq. (44), but where now all fields are subject to the periodic 
boundary conditions (139). 

Z(B) = N\detW(B)); e l J 2 [det(-D,(B)D ll (B + a))] a= s 



S.T 
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X 1 + 



E 

1=1 



s 



in1 



I \ 



l\ 



> exp (~JGj) 

V 2 / j =D 



(142) 



We can consider strong background fields, so that gB^ is large (as compared to 
Aq CD ), while a s = g 2 /An in that strong background is small at all distances. 
Moreover, it was shown that a s is frozen at large distances [17]. In this case 
Eq. (142) is a perturbative sum in powers of g n , arising from the expansion in 

(ga») n . 

In what follows we shall discuss the Feynman graphs for the free energy F(T), 
connected to Z(B) via 



F{T) = -T\n{Z{B)) B . 



(143) 



As will be seen, the lowest order graphs already contain a nontrivial dynamical 
mechanism for the deconfmement transition, and those will be considered in 
the next subsection. 

6.2 The lowest order gluon contribution 



To lowest order in ga^ (keeping all dependence on gB^ explicit) we have 



Z = e~W/ T = N'(exp(-F (B)/T)) B , 



(144) 



where using Eq. (142) Fq(B) can be written as 



^F (B) = i In det^r 1 - \ndet(-D 2 (B)) 

J. Zi 



foo OO 



(145) 
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In Eq. (145) Sp implies summing over all variables (Lorentz and color indices 
and coordinates) and is a regularization factor as in Eq. (134). Graphically, 
the first term on the r.h.s. of Eq. (145) is a gluon loop in the background field, 
while the second term is a ghost loop. 

Let us turn now to the averaging procedure in Eq. (144). With the notation 
<f = —F (B)/T, we can exploit in Eq. (144) the cluster expansion[14] 

(exp</?) B =exp (^T((<p n }}^J 

= ex P {( V } B + - (^) B ] + 0(y, 3 )}. (146) 

To get a closer look at (ip) b we first should discuss the thermal propagators of 
the gluon and ghost in the background field. We start with the thermal ghost 
propagator and write the FSR for it [11] 

oo oo 

(-DX = {x\ J dte tD2 W\y) = J dt(Dzr xy e- K $(x,y). (147) 



Here $ is the parallel transporter in the adjoint representation along the tra- 
jectory of the ghost: 

$(x,y) = Pexp(ig J B^(z)dz^) (148) 

and (Dz)™ is a path integration with boundary conditions imbedded (denoted 
by the subscript (xy)) and with all possible windings in the Euclidean temporal 
direction (denoted by the superscript w). We can write it explicitly as 



(Dz)™ = hm [ sv ' 

V N ^rntl (4^) 2 

x d 4 p 



N 



ip C( m ) -( x ~v)- n /%4 

\rn=l / 



(149) 
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Here, ((k) = z(k) — z(k — 1), Ne = t. We can readily verify that in the free 
case, Bfj, = 0, Eq. (147) reduces to the well-known form of the free propagator 



f ( 2 (m) 



-v d 4 p 



(2k) 



x exp \ip (J2 C(m) - (x - y) - n[35^ 4 ) 

oo 

E / exp [—p 2 t — ip(x — z) — ip 4 n[3 



d 4 P 

(2k) a 



(150) 



with 



d((m) = 



d((m) 
(Aire) 2 



Using the Poisson summation formula 
1 



2tt 



E exp(ip 4 n(3) = E 5(p A f3-2nk) 



(151) 



n=0,±l,±2.. 



fc=0,±l, 



we finally obtain the standard form 

Td 3 p exp[— ipi(x — y)i — i2nkT(x4 — y^)] 



ixy 



E 



fc=0,±l,.. 



(2k)* 



pj + (2nkTf 



(152) 



Note that, as expected, the propagators (147) and (152) correspond to a sum 
of ghost paths with all possible windings around the torus. The momentum 
integration in Eq. (149) asserts that the sum of all infinitesimal "walks" ((m) 
should be equal to the distance (x — y) modulo iV windings in the compacti- 
fied fourth coordinate. For the gluon propagator in the background gauge we 
obtain similarly to Eq. (147) 



G xy = I dt(Dzr xy e- K <S> F (x,y), 



(153) 



where 



$ F (x,y) = P F Pexp -2ig J F(z(r))dr J exp iig J B^dz^ 



(154) 
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The operators P F P are used to order insertions of F on the trajectory of the 
gluon. 



Now we come back to the first term in Eq. (146), (</?)#, which can be repre- 
sentated with the help of Eqs. (147) and (153) as 

'1 



(<p)b = J j((t)d*x(Dzr xx e^ 



^tr(^ F (x,x)) B - (tr®(x,x)) B 



, (155) 



where tr implies summation over Lorentz and color indices. We can easily 
show [11] that Eq. (155) yields for £> M = the usual result for the free gluon 
gas: 

F (B = 0) = -Tip(B = 0) = -(N 2 - m^f- (156) 
6.3 The lowest order quark contribution 

Integrating over the quark fields in Eq. (140) leads to the following additional 
factor in Eq. (142) 

det(m + D(B + a)) = [det(m 2 - D 2 {B + a))} 1 ' 2 . (157) 

In the lowest approximation, we may omit a M in Eq. (157). As a result we get 
a contribution from the quark fields to the free energy 

oo 

= -ilndet(m 2 - D 2 {B)) = - l -Sp J ^e^^^, (158) 

o 

where Sp has the same meaning as in Eq. (145) and 

D 2 = (D^) 2 = Dl(B) - gF^a^ = D 2 — gaF; 

<V = + 4 {l^lv - lulu) • (159) 
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Our aim now is to exploit the FSR to represent Eq. (158) in a form of the 
path integral, as was done for gluons in Eq. (147). The equivalent form for 
quarks must implement the antisymmetric boundary conditions pertinent to 
fermions. We find 



1 1 °° 



(160) 



where 



and 



/ 



W ff (C n ) = P F P A exp 



ig J A^dz^ exp g (crF) , 



(Dz\ 



xy 



\ (4tt£) 



n=0,±l,±2,. 



cxp 



/ N \ 

ip CM ~(x-y)- 

\m=l / 



.(161) 



It can readily be checked that in the case = the well known expression 
for the free quark gas is recovered, i.e. 



Fq (free quark) 



III 
180 



N c V z T A -n fl 



(162) 



where n/ is the number of flavors. The derivation of Eq. (162) starting from 
the path-integral form (160) is done similarly to the gluon case given in the 
Appendix of the last reference in Ref. [11]. 

The loop C n in Eq. (160) corresponds to n windings in the fourth direction. 
Above the deconfinement transition temperature T c one sees in Eq. (160) the 
appearance of the factor 



Q = P exp 



ig J B i {z)dz i 
o 



(163) 
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For the constant field B 4 and Bi — 0, i — 1, 2, 3, we obtain 

W = -^E^-(-ir. (164) 



This result coincides with the one obtained in the literature [36]. 



7 Discussion and conclusions 



Three basic approaches to QCD which are largely used till now are: i) lat- 
tice simulations ii) standard perturbation theory, and iii) OPE and QCD 
sum rules[41]. The two latter methods are analytic and have given enormous 
amount of theoretical information about the high-energy domain, where per- 
turbative methods are applicable, and about nonperturbative effects both in 
the high and low energy regions. 

These methods have their own limitations. In particular, the standard pertur- 
bation theory is plagued by the Landau ghost pole and IR renormalons and 
slow convergence, which necessitates the introduction of methods, where sum- 
mation of perturbative subseries can be done automatically and the Landau 
ghost pole is absent. The QCD sum rules are limited by the use of only a few 
OPE terms, while the OPE series is known to be badly convergent and at best 
asymptotic. 

One of the great challenges of QCD is to have a tractable analytic treatment of 
it. In particular, the improvement of the standard perturbation theory and the 
search for a systematic approach to nonperturbative phenomena are impor- 
tant objectives. The methods presented here in the present paper, commonly 
entitled The Fock-Feynman-Schwinger Representation, are meant to exactly 
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do this. The main advantage of the FSR is that it allows to treat both per- 
turbative and nonperturbative configurations of the gluonic fields. 

In case of purely perturbative fields the FSR yields a simple method of summa- 
tion and exponentiation of perturbative diagrams[15]. Nonperturbative fields 
are introduced in the FSR naturally via the Field Correlator Method[29,56]. A 
recent discovery on the lattice of the Gaussian correlator dominance (Gaussian 
Stochastic Model) (see Ref. [46] for discussion and further references.) makes 
this method accurate (up to a few percent). There is another very important 
result of taking into account nonperturbative fields in the QCD vacuum: this 
fact allows to develop perturbation theory in the nonperturbative background 
- which is realistic unlike the standard perturbation theory. It contains no 
Landau ghost poles and IR renormalons[17]. 

As two applications of the FSR we have considered the problem of collinear 
singularities and finite temperature QCD. As an important special feature we 
should stress the absence in the background perturbation theory of all IR and 
collinear singularities pertinent to standard perturbation theory. This feature 
discussed here, opens new perspectives to the application of FSR to high- 
energy QCD processes. Although not discussed here, the FSR can readily be 
extended to treat deep inelastic scattering, Drell-Yan and other processes, in- 
cluding the fundamental problem of the connection between constituent quark- 
gluon model and parton model. 

Moreover, as shown here the FSR can be used to describe QCD at nonzero tem- 
perature at above and around phase transition point. As was shown before [11- 
13] the vacuum is predominantly magnetic and nonperturbative above T c . 
Therefore methods based on FSR are working well in this region. In conclu- 
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sion, since the FSR is replacing field degrees of freedom by corresponding 
quantum mechanical ones it has the advantage, that the results can often be 
interpreted in a simple and transparant way. It has been applied with success 
to both Abelian and non-Abelian situations. We have found that the FSR is 
a powerful approach for studying problems in QCD. 
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